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Abstract. We calculate stochastic gravitational wave 
background produced by extragalactic merging binary 
white dwarfs at the LISA frequencies 10 -3 — 10 -2 Hz with 
account of a strong evolution of global star formation rate 
in the Universe recently established observationally. We 
show that for the observed global star formation history 
and modern cosmological models the extragalactic back- 
ground is an order of magnitude smaller than the mean 
Galactic value. An early star formation burst at high red- 
shifts can bring it at a higher level but still a few times 
lower than the mean Galactic one. 
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1. Introduction 

Gravitational waves (GW) from unresolved binary stars 
in the Galaxy forms a confusion limit in the frequency 
domain of LISA laser space interferometer (Bender et al. 
1996) 1CT 4 - 1CT 1 Hz that effectively adds to the detec- 
tor's noise. Different aspects of this background have been 
studied earlier by Rosi and Zimmerman (1976), Lipunov 
& Postnov (1987), Lipunov, Postnov & Prokhorov (1987), 
Hils, Bender & Webbink (1990), Bender & Hils (1997), Gi- 
ampieri & Polnarev (1997), Postnov & Prokhorov (1998). 

The importance of the knowledge of binary confusion 
limit within the LISA frequency band is dictated by many 
factors - for example, it limits the angular resolution of 
the detector (Cutler 1997) and restricts possibilities of ob- 
serving possible relic GW backgrounds (see, e.g., discus- 
sion in Grishchuk 1997). Since LISA with its three arms 
is expected to operate as essentially single interferome- 
ter (Schutz 1997), it can only see noise sources above its 
intrinsic noise h rms . In our previous paper (Postnov & 
Prokhorov 1998) we focused on the stochastic GW noise 
from galactic binaries and showed that, depending on (un- 
known) galactic binary white dwarf merger rate, this noise 



becomes lower than the expected LISA rms noise at fre- 
quencies / ~ 3 x 1CP 2 Hz B We concluded that to be 
detectable by LISA, relic GW backgrounds should be as 
high as SlGwh\ m > 1CT 8 at 10~ 2 Hz. 

Clearly, the galactic binary GW noise, tracing the 
distribution of stars in the Galaxy, should be highly 
anisotropic (of order of one magnitude higher in the di- 
rection of the Galactic center; see calculations by Lipunov 
et al. 1995) and thus modulated with LISA turning in its 
orbit. This distincitve feature of the galactic stochastic 
background can in fact be used to measure it by space- 
born interferometers (Giampieri & Polnarev 1997). This 
is not the case for a GW-background produced by ex- 
tragalactic binaries, which becomes important at frequen- 
cies where galactic binaries do not form continuous noise 
(/ > 3 x 10 -3 Hz for 1-year integration). 

A crude estimate shows that the isotropic extragalac- 
tic background is expected to be an order of magnitude 
smaller than the average GW noise from galactic binaries 
(see Lipunov et al. (1987, 1995), Hils et al. (1990)). All 
estimates of this background made so far, however, have 
not taken into account the fact of a strong global star for- 
mation evolution recently revealed by different astronomi- 
cal observations (see the results of Canada-France-Hawaii 
survey of far galaxies with z < 1 (Lilly et al. 1996), the 
analysis of galaxies with Lyman jump at z > 2 (Madau et 
al. 1996), and the analysis of the HST deep field survey of 
galaxies with 1 < z < 2 (Connoly et al. 1997)). According 
to these studies, the global star formation rate strongly in- 
creases with redshift, peacking at z ~ 2 (see Fig. 1 below 
reproduced from Connoly et al. (1997), with points from 
Madau et al. (1996) corrected for extinction as in Boyle 
& Terlevich (1997)). 



1 This limiting frequency is derived assuming integration 
time T and frequency / such that fT = 1; for a fixed inte- 
gration time (say, 1 year) the frequency resolution is A/ ~ 
1/T ~ 10~ 7 Hz. Correspondingly, the limiting frequency for 
the galactic white dwarf binary confusion limit would decrease 
as (fT)~ 3 / 8 , typically by 50 times at 1 mHz for 1 year integra- 
tion time. 
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In a Euclidean space, the remote sources would con- 
tribute more and more to the energy density with distance 
r, as oc r, leading to the photometric paradox; in a FRW 
expanding Universe this is not the case. However, the ef- 
fects of proper density evolution of galaxies (star forma- 
tion) with redshift, together with the proper evolution of 
the sources inside galaxies, could compensate for the cos- 
mological energy density dilation. So a priori it is not at all 
clear if the remote extragalactic binaries can be neglected 
in studying extragalactic binary GW background. This is 
the motivation to our calculations. 

2. GW noise from sources with changing frequency 

We start with a short description of how a GW back- 
ground is formed by a population of frequency-unresolved 
sources with changing frequency. The amplitude of the 
stochastic GW background is conventionally measured in 
terms of the energy density per logarithmic frequency in- 
terval related to the critical density to close the Universe 
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where p cr = 3Hq/8it w 1.9 x 10~ 29 /i 2 00 g cm~ 3 , h wo = 
-f/o/100 [km/s/Mpc] is the present value of the Hubble 
constant; we shall use geometrical units G = c = 1 
throughout the paper. As in the case of deterministic GW 
signals, one may equally use the characteristic dimension- 
less amplitude of metric variations h c , which determines 
the signal-to-noise ratio on the detector (S/N) — h c /h rms . 
This amplitude is connected with Slew through h c (f) = 

(l/27T)(ff //)n^. 

To calculate GW noise produced by some frequency- 
unresolved sources, we need know the number of sources 
per logarithmic frequency interval. At the first glance, this 
would require knowing the precise formation and evolution 
of sources. However, when only GW carries away the angu- 
lar momentum from the emitting source (as is the case for 
coalescing binary white dwarf forming the high-frequency 
part of the galactic binary GW background) , the problem 
becomes very simple and physically clear. 

Energy loss causes the change in the GW frequency 
u> = 2nf of emitting objects. In the case of a bi- 
nary star with masses of the components M\ and Mi 
in a circular orbit of radius a the negative orbital en- 
ergy, E orb = -M x M 2 /2a ~ MflU (here M = (M 1 + 
M 2 ) 2/5 (MiM 2 ) 3/5 is the so-called "chirp mass" of the sys- 
tem) is being lost, and the orbital frequency (hence, the 
GW frequency, which is twice the orbital one) increases. 

To a good approximation, the conditions of star forma- 
tion and evolution of astrophysical objects in our Galaxy 
may be viewed as stationary. This is true at least for last 
5 billion years. Let the formation rate of GW sources be 
1Z. Inside the LISA frequency range, 10~ 4 — 1CP 1 Hz, only 
coalescing binary white dwarfs (WD) and binary neutron 



stars contribute. Even if binary neutron stars coalesce at a 
rate of 1/10000 yr in the Galaxy, their total number at any 
moment of time still should be much smaller than of the 
WD binaries, and we restrict ourselves to considering only 
binary WD. As we showed in our previous paper (Post- 
nov & Prokhorov 1998), the galactic rate of binary WD 
mergers may be restricted by the statistics of SN type la 
explosions, which are much more numerous than the ob- 
served close binary WDs. However, as we show below, the 
ratio of the extragalactic to galactic GW binary noise is 
almost independent of the exact value of TZ. In real galax- 
ies 1Z may be the function of time, which will be accounted 
for by the star formation history SFR(i). 

At some moment of time the number of sources per 
unit logarithmic frequency interval is 



dN(t)/dlnf = N(f,t)=Tl(t)(f/f). 



(2) 



so the total energy emitted in GW per second per unit 
logarithmic frequency interval at frequency / by all such 
sources in the galaxy is 



(!) dE GW /(dtd\nf) 



L(f) GW = L(f) GW N(f,t) 
L(f) GW K(t) x (///) , 



(3) 



where L(f ) G w is the GW luminosity (erg/s) of the typ- 
ical source at frequency /, L(f) G w ^ f 10 ^ 3 f° r binary 
systems. 

For an isotropic background we find 



tt G w(f)Pcr = L(f) GW /(4TT(r) 2 ) 



(4) 



where (r) is the inverse-square average distance to the 
typical source. Strictly speaking, this distance (as well as 
the binary chirp mass M) may be a function of frequency 
since the binaries characterized by different M may be 
differently distributed in the Galaxy. Since the real dis- 
tribution of binaries of different mass in the Galaxy is 
poorly known, we can take, as an estimate, the mean pho- 
tometric distance for a spheroidal distribution in the form 
dN oc exp[— r/ro] exp[— (z/zq) 2 ] (r is the radial distance 
to the Galactic center and z the hight above the galactic 
plane) with r$ — 5 kpc and zq — 4.2 kpc with (r) « 7.89 
kpc. 

For binary stars the energy reservoir radiated in 
GW is their orbital binding energy, which is a power 
law of the orbital frequency: E oc / 2//3 . Hence the 
frequecy change (///) may be found from the equa- 
tion dE/dt = (2/3)E(f I f) By energy conservation law 
dE/dt = (dE/dt) G w + (dE/dt)... where (. . .) means other 
possible loss/gain of energy. Finally, we obtain 



(///) = (2/3)E/((dE/dt) GW + (dE/dt)...) 
and 

L(f) GW = {2/3)KE- 



1 + 



(dE/dt)... 
(dE/dt) G w 



(5) 



(6) 
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Remarkably, if GW is the dominant source of energy re- 
moval, the resulting GW stochastic background depends 
only on the source formation rate: 



tt G w(f)Pcr = (2/3)KE/(4ir(r) 2 



(7) 



For the frequency range under consideration the as- 
sumption that only GW emission drives the orbital evolu- 
tion of a close binary WD is a perfect approximation. The 
same relates to the assumption tacitly made that the orbit 
of all such binaries is circular. As we are interested in the 
blue end of this frequency range, we need not to consider 
much more complicated evolution of cataclysmic binaries 
whose orbital periods are always higher than a few tens of 
minutes. 

3. GW noise from unresolved extragalactic binary 
stars at LISA frequencies 

Let the detector's frequency be /. Then sources in galaxies 
within the comoving volume element dV(z) at redshift 
z produce an isotropic GW background with the energy 
densidy near the detector (at z = 0) 



p cr dfl GW (f) 



( HP 
\4irD(z) 2 



n{z)dV{z) , 



(8) 



where L(f) is the GW luminosity of one galaxy at the 
comoving frequency /' = f(l + z), D(z) is the photo- 
metric distance, n{z) = n G F(z)(l + z) 3 is the comoving 
density of galaxies at redshift z with factor F(z) describ- 
ing the galactic density evolution, and n G is the present 
day density of galaxies. In the case of no proper evolution 

F{z) = 1 \. 

Transiting in a standard way from integration over 
proper volume to integration over cosmic time (e.g. Wein- 
berg 1972) we arrive at 



p cr fl G w(z*) 



n G I F{z) 
h 



£(/(! + *)) 
1 + z 



dt , 
— — dz . 

dz 



(9) 



where z* marks the beginning of the star formation in the 
Universe. Noticing that L(f') oc E or b{f(l + zj) scales with 
redshift as (1 + z) 2 / 3 , using Eq. (0), and assuming that 
binary WD merging rate 1Z per proper volume changes 
with time only due to galactic density evolution F{z), we 
finally obtain the ratio of extragalactic to galactic GW 
energy density produced by coalescing binary WD: 



^G\v(z*) 



GW 



47r(r 2 )riG 



F(z) dt 
o {l + z)-^ 3 Tz 



dz 



(10) 



The present-day galactic density can be rewritten in 
terms of fib, the baryon fraction in stars relative to the 
critical density, assuming the stellar mass of the typical 
galaxy to be M G = 10 n M Q : 

fi b 



n G « 0.013(Mpc~ 



0.005 



noo • 



(11) 



Here fib is normalized to the observed present-day baryon 
density of luminous matter (fi(, = 0.003 inside bulges of 
spiral and in elliptical galaxies, fib = 0.0015 in disks of 
spiral galaxies; see Fukugita et al 1996, 1997 for further 
details). For hioo — 0.7 the numerical coefficient in (11) 



is in a good agreement with the measured local galactic 
density 0.0048Mpc" 3 (Loveday et al. 1992). 

The time-redshift relation in a FRW cosmological 
model with arbitrary matter (Qm) an d vacuum (fl\) den- 
sity is 



dt = 



1 



Ho (1 + z)y/(l + z) 2 (l + fl M z) - z(2 + z)fl A 



(12) 



(notice that flM includes both visual baryon density fib 
and dark matter). 

For a flat (CIm = 1) FRW universe without cosmologi- 
cal term (Qa = 0) dt/dz = (l + z)~ 5 / 2 and without source 
evolution (F(z) = lV,the integrand in Eq. ( |lO| ) turns into 



(1+*) 



-17/6 



so 



thatH 



£Igw{z*) 



GW 



0.03 



fin 



0.005 



tioo 



<r) 



lOkpc 



[i-(i + zT 11/6 ] 



(13) 



Clearly, the extragalactic background is small in compar- 
ison with the galactic one if no cosmological evolution of 
global star formation rate is included. 

Now we wish to take into account the evolution of WD 
merging rate with time (redshift). Let G(t) be the time 
dependence of the event rate of interest after a (5-function- 
like star fromation burst. This function can be considered 
as a Green function for arbitrary law of star formation 
SFR(t). Then the comoving event rate at any time can be 
calculated as 



K G (t) 



X 

J SFR(r)G(t - r)dr 



(14) 



*(*.) 



(here t(z*) is the initial moment of star formation). The 
event rate per proper time observed from a layer dz at 
redshift z is 



K(z) = n{z)K G {z)dV{z). 



(15) 



Clearly, the evolution of the comoving galactic density 
n(t(z)) can be considered in the same way, but since we de- 
rive from observations the comoving luminosity density as 
a function of redshift, it is impossible to separate the star 
formation rate inside galaxies and the comoving galactic 
density evolution without making special model assump- 
tions. We will assume that expressing the comoving den- 
sity of galaxies through the luminous baryon density fib 



2 Eq. (14) in paper (Postnov & Prokhorov (1998)) was cal- 
culated for F(z) — (1 + z) :i , so it should be substituted by this 
equation for non-evolving sources. 
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( |ll| ) reflects real global SFR evolution without knowing 
how precisely has the comoving density of different types 
of galaxies changed with time. This remains valid until Qb 
is assumed constant. For example, fif, may change in mod- 
els with non-zero cosmological constant at stages when its 
influence on the expansion of the Universe becomes dom- 
inant. 



R = 1 + 12.2z + -3.65z 3 + 0.291z 5 




Fig. 1. The evolution of global star formation rate in the Uni- 
verse from Connolly et al (1997) normalized such that the 
present star formation rate in solar masses per year per a 
10 11 M o galaxy is unity. Open squares at large redshifts repre- 
sent data from Madau et al. (1996) corrected for extinction as 
in Boyle and Terlevich (1997). 

The observed star formation rate evolution SFR(z) 
were taken from Connolly et al. (1997). In fact, the ab- 
solute value of star formation rate is deduced from the 
observed UV luminosity density, which requires knowledge 
of the initial mass function of stars and models of stellar 
evolution (Madau et al. 1997). For the relative star forma- 
tion these requirements become weaker (at least so far as 
we assume the same parameters of the star formation and 
evolution at any time). We reproduce this function (nor- 
malized such that SFR(z = 0) = 1) in Fig. |l] with data 
from Madau et al. (1996) (open squares) being corrected 
for extinction as in Boyle and Terlevich (1997) (then the 
global SFR coincides with the evolution of the comoving 
density of quasars). For numerical calculations this SFR 
can be approximated by the polynomial 

SFR(z) m 1 + 12.22 - 3.65z 2 + 0.291z 3 . (16) 

Fig. ^| shows the dependence of the binary white dwarf 
coalescence rate per a IO^Mq galaxy on redshift calcu- 
lated using Eq. ( |l4| ) in units of the (precisely unknown) 
galactic rate lZ w d for different representative cosmologi- 
cal models: the standard (SCDM: il tot — 1,0a = 0), 
a flat with the cosmological constant (AC DAI : S! fot = 
l,f2A = 0.7, and an open without the cosmological con- 
stant (OCDM: Qtot = 0.2, A = 0). The Green function 
G(t) for binary WD coalescence rate was taken from our 
numerical calculations (Lipunov & Postnov 1988; Lipunov 
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Fig. 2. The dependence of binary WD coalescence rate 
(normalized to the present Galactic rate) on redshift z for 
three representative cosmological models: a standard one 
(yitot = 1,^a = 0), a flat with cosmological constant 
(Qtot = 1,^a = 0.7, and an open without cosmological con- 
stant (n to t = 0.2, fi A = 0). 

et al. 1995; Lipunov, Postnov & Prokhorov 1996) which 
to a good approximation is 

G W = TT^P' (17) 

where tu is the Hubble time, k is a numerical fac- 
tor, R is a constant calculated assuming that the present 
Galactic binary WD coalescence rate is the mean value for 
a constant star formation: 

iff 

n = l^U-) I G{T)dT (18) 

i(*„) 

It is seen from this figure that the choice of the cos- 
mological model affects this function insignificantly. 

After having obtained TZg(z) 1 we can easily calcu- 
late, substituting F(z) — TZq( z ) into Eq. JlC|), the ratio 
Qgw(< z„) /Q,qw(z — 0) for different cosmological mod- 
els. The results are presented in Fig. ||for the same models 
as in the previous Figure. For comparison, the thick solid 
line shows the value Q,gw(< z*)/flGW assuming a con- 
stant star formation rate in the Universe (Eq. fl3|). 

Our consideration of the stochastic GW background 
formation in the process of extragalactic binary WD merg- 
ing remains valid until redshifts where " boundary effects" 
of frequency distribution of binary stars are unimportant. 
In the frequency range f b s = 10 -3 — 1CU 2 Hz this is al- 
ways the case because the maximum proper frequency of 
the noise f em = / fcs(l + z) < 10/ em < 0.1 Hz for z* < 6. 
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Fig. 3. The ratio of the energy density of the extragalactic 
stochastic GW background formed by merging binary WD 
Qgw(< z*) to the mean Galactic level at frequency / = 0.01 
Hz as a function of parameter z* , the epoque of the inital star 
formation in the Universe. The global star formation rate his- 
tory has the form jl^). Cosmological models are the same as 
in Fig. |. 



At present, the lack of observational data on the SFR 
behaviour at high redshift does not allow us to make more 
robust estimates. We conclude that unless the global star 
formation rate continues increasing with redshift at z > 3, 
the extragalactic GW background energy density ^IqW ~ 
is ten times smaller than the mean galactic value in 
the LISA frequency range 10~ 4 — 10 _1 Hz. 

The galactic binary GW noise will be modulated by 
the LISA orbital motion, whereas the extragalactic one 
will not. If the latter is comparable with galactic values at 
some galactic latitudes, it can impede detection by LISA 
of some interesting relic cosmological GW backgrounds 
(e.g. Grishchuk 1997), and specific statistical features of 
the relic GW should be used to separate them against the 
noise from astrophysical sources. 

The authors acknowledge D. Hils and P.L. Bender for 
correcting an error in our calculations. The work was 
partially supported by Russian Fund for Basic Research 
through Grant No 98-02-16801, the INTAS Grant No 93- 
3364, and by NTP "Astronomy" (project 1.4.4.1) of the 
Ministry of Science of Russia. 



4. Conclusions 



References 
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background, together with independent studies of far-IR 
background, would help revealing the true star formation 
history at high redshifts. 
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